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Abstract 

Light localization is controlled at a scale of A/ 10 in the harmonic regime from the far field domain 
in a plasmonic nanoantenna. The nanoantenna under study consists of 3 aligned spheres 50 nm in 
diameter separated by a distance of 5 nm. By simply tuning the orientation of an incident plane 
wave, symmetric and antisymmetric mode-balancing induces a strong enhancement of the near 
field intensity in one cavity while nullifying the light intensity in the other cavity. Furthermore, it 
is demonstrated that the dipolar moment of a plasmonic particle can be fully extinguished when 
strongly coupled with a dimer of identical nanoparticles. Consequently, optical transparency can 
be achieved in an ultra-compact symmetric metallic structure. 

PACS numbers: Valid PACS appear here 
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Metallic nanostructures offer the opportunity to strongly focus light and to enhance light- 
matter interactions at the nanometer scale via the excitation of localized surface plasmon 
resonances. In 2002, Stockman et al. demonstrated the dynamic control of the field enhance- 
ment at nanometer scale in arbitrary nanostructures pQ by adjusting the phase and the po- 
larization of an excitation pulse. Active control of field localization has been experimentally 
demonstrated through the use of an iterative learning algorithm to shape the femtosecond 
excitation pulse [2j. Deterministic structures such as arrays of spherical nanoparticles have 
been investigated by Koenderinck et al. [3j to serve as unique and reproducible lithographic 
masks for imprinting different near field patterns. They varied the angle of incidence and 
polarization of the continuous excitation beam in order to produce a variety of patterns in 
a photosensitive substrate. These studies convincingly demonstrated that incident beam 
shaping can induce different field distributions in a single metallic nanostructure. 

In 2006, Le Perchec et al. controlled the localization of light in coupled slits milled in a 
metallic substrate [4] . The slits had a thickness of 200 nm and were separated by a distance 
of 500 nm. When illuminated in oblique incidence in the infrared spectrum, they coupled 
incident propagating light to both symmetric and anti-symmetric modes. By the use of a 
second identical frequency propagating wave in opposite incidence, they could switch the 
near field enhancement from one slit to the other via phase adjustments. This simple system 
provided a means to conceive subwavelength optical switch components controlled from the 
far field region. Very recently, Volpe et al. [5] investigated spatial phase modulations of high 
order beams such as Laguerre-Gaussian beams in order to control the near field in plasmonic 
nanostructures. They considered two 50 nm gaps formed by three 500 nm aligned gold bars 
and demonstrated that spatial phase modulation of the incident beam enables switching light 
on and off in either of the two sites. The aim of this paper is to theoretically study a linear 
trimer of identical spherical nanoparticles to tune the localization of light at the nanometer 
scale (^55 nm). In particular, we investigate the near field enhancement occuring between 
gold nanospheres 50 nm in diameter. The three spherical nanoparticles are arranged to 
form two identical cavities in which the near field can be strongly enhanced. This provides 
a simple system with two coupled nanocavities separated by a distance of about a tenth of 
the incident wavelength. It will be shown that this a priori symmetric system enables the 
concentration of light in a single cavity while employing a single illumination plane wave. 

We investigate a system consisting of three identical gold nanospheres placed along the 
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FIG. 1: Layout of the numerical experiment. Three identical gold nanospheres 50 nm in diameter, 
embedded in a dielectric background of refractive index no = 1.5, are placed along the z-axis. The 
system is illuminated by a plane wave with a TM polarization with an angle of incidence 9 between 
the propagation direction and the z-axis. Calculations are performed using generalized Mie theory 
combined with an analytic multiple scattering formulation. 

x-axis with a diameter of 50 nm and embedded in a polymer substrate with a refractive 
index of 1.5 (Fig. [I]). The choice of the structure is motivated by the study of Bidault et 
al. [6] where well-defined groupings of Au nanoparticles with controlled nanometer spacing 
were synthesized by hybridizing monoconjugated DNA in a polymeric material of refractive 
index n = 1.5. This bottom-up technology opens up a way to conceive coupled cavities 
made of nanoparticles separated by a few nanometers, currently unobtainable via classical 
lithography techniques. 

This model system is illuminated by a plane wave with TM polarization (the magnetic 
field is perpendicular to the plane presented in Fig. [j]) and an angle of incidence, 0, tunable 
between (propagation along the z-axis) and 90 degrees. The strong contrasts of near field 
intensities at scales much smaller than the incident wavelength require a full electromagnetic 
study based on solving the Maxwell equations. Rigorous Lorentz-Mie theory is employed in 
order to accurately calculate the near field enhancement in the vicinity of a single spheri- 
cal metallic nanoparticle. With this technique, electromagnetic fields inside the sphere and 
in the surrounding medium are expanded using basis sets of Vector Spherical Harmonics 
(VSHs), which allow an analytic satisfaction of the boundaries conditions. The expansion 
of an arbitrary vector field is given in Eq. ([!]) in terms of the Vector Spherical Harmonics. 
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The VSHs provide an angular description of the vectorial field and are defined in terms of 

) . The En,m ' are their respective coef- 
ficients that describe the radial dependency of the field according to the medium properties. 
Interactions between particles are taken into account by invoking the translation/addition 
theorem in order to solve the Foldy-Lax multiple scattering equations [8]. The truncation 
of the basis sets, necessary for a numerical solving of this system of equations, is deter- 
mined from the size parameter of the spheres and the strength of electromagnetic couplings 
between the particles, and results in a finite dimensional T-matrix describing the electro- 
magnetic response of the entire system [9] . In this study, high multipolar orders are required 
to accurately reconstruct the field which will be taken into account by imposing n max = 10 
in Eq. ([T]). The dispersion of the dielectric constant of gold is interpolated from the Palik 
database HOI. 
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We first investigate the far field response of the trimer by rigorously calculating the spec- 
tral response of the extinction cross-section for different angles of incidence (Fig. [2)2). For an 
incident angle 6 = 0° (blue line) corresponding to an illumination along the z-axis, only one 
maximum appears at the wavelength of 560 nm, which is almost the wavelength of resonance 
of the isolated sphere (red line). Assuming an incident electric field along the x-axis, it can 
be deduced that this resonance is due to in-phase plasmon oscillations occuring transver- 
sally to the chain axis resulting in weak couplings between the individual nanoparticles 
resonances. At normal incidence (6 = 90°: black line), a stronger resonance occurs at the 
wavelength of 665 nm, corresponding to the well-known red-shifted coupled resonance of in- 
phase longitudinal modes [llj. Fig.[2|) illustrates the strong near field enhancement between 
the spheres when the trimer is illuminated at A = 665 nm. When the chain is illuminated 



at normal incidence, only in-phase longitudinal oscillations of the polarization moments are 
excited. These oscillations induce surface charges of opposite sign in the dielectric inter- 
particle gap responsible for the strong electric field enhancement. Such longitudinal modes 
have been studied extensively in the context of dimers of metallic spheres and nanoshells 
[TTl - fT6] . We remark that for the case of normal incidence, there exists a weak second peak 
at higher frequency which is made possible thanks to quadrupolar response contributions. 
At intermediate angles of incidence (like the case 9 = 45° displayed in Fig. [2)2), we observe 
two peaks at frequencies close to those observed respectively at 9 = 0° and 9 = 90°. 
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FIG. 2: (a) Extinction cross-section of a chain of three gold particles of diameter D = 50 nm (red 
line) separated by w = 5 nm in the full electromagnetic approach. Different angles of incidence 
are displayed: = 0° (blue line), 9 = 45° (green line) and 9 = 90° (black line). The isolated sphere 
extinction cross-section is also displayed (red line), (b) Electric field intensity ||£'|| 2 /||£'^ ric || 2 map 
in log scale for a trimer illuminated at A = 665 nm under normal incidence 

Table. [T] illustrates the dipolar eigenmodes that may be excited in a linear trimer. Al- 
though the dipolar approximation fails to describe accurately the resonance when the parti- 
cles are strongly coupled, it enables a clear distinction between modes that are perpendicular 
(transverse) and parallel (longitudinal) to the chain axis [17J. From the extinction cross- 
section displayed in Fig. [2)2, one can observe that only the in-phase modes denoted Tl and 
LI are clearly visible. This observation may seem surprising since except for illumination 
at 9 = 0° or 9 = 90° where the symmetries of the configuration impose in-phase modes, 
coupling between the incident field and opposite-phase modes may be possible when illumi- 
nating the system at oblique incidence. Inspection of the opposite-phase modes in Table. [T] 
shows however that the total dipolar moment involved in these modes is weaker than that 
occurring for in-phase modes. Consequently, these modes are not easily coupled to the far 
field and cannot be observed in extinction cross-section plots. 
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TABLE I: Possible dipolar eigen modes for a linear chain of three plasmonic particles 



Transverse eigenmodes 


Tl | t t 


In-phase mode 


T2 i t t 
T3 t I t 


Opposite-phase 
modes 


Longitudinal eigenmodes 


LI -> -> -> 


In-phase mode 


L2 < > -> 

L3 < > 


Opposite-phase 
modes 



Let us emphasize that we are interested in the vanishing of the near field enhancements 
in one of the cavities. This purpose requires an excitation via the far field of longitudinal 
opposite-phase modes, modes that are associated with a weak dipolar moment. Fig. [3)2 
displays the near field intensity enhancement calculated at the center of the cavities as a 
function of the incident wavelength, when incident angle 9 is equal to 45 degrees. These plots 
show that the field intensities in the two cavities are quite different, but most interestingly 
they reveal the counterintiutive occurence that at A = 595 nm (green line) a rather strong 
field intensity enhancement in gap 2 is concurrent with a near vanishing of the field in gap 
1 (ratio larger than 500). This phenomenon is even more conspicuous in the near field 
intensity map displayed in Fig. [3j> in logarithmic scale that clearly shows that at A = 595 
nm, the light intensity is confined in only one site. This is the most striking result of this 
study: the sub-wavelength cavity 1 is turned-off and the light intensity is almost entirely 
localized inside cavity 2 while both cavities are separated by only 55 nm. It appears that 
this phenomenon cannot be observed from the far field response since the extinction cross 
section (green line in Fig. [2)2) does not exhibit any notable feature. We next move towards 
the full understanding of the longitudinal eigen modes involved in this phenomenon, which 
will require the extraction of longitudinal terms in the multipole formalism. 

The extraction of the longitudinal components of the field requires that we focus attention 
on the Scalar Spherical Harmonics that are defined in Eq. ^ in terms of the associated 
Legendre polynomials P™(cos#). 

Y nm (9, 0) = c n , m P™ (cos 9) exp(zm0) (3) 
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FIG. 3: (a) Intensity enhancement in cavity 1 (red line) and cavity 2 (black line) for 9 = 45°. The 
vertical line indicates the wavelength of 595 nm where the intensity contrast is maximum, (b) 



Electric field intensity, \\E\ 



map in log scale at A = 595 nm and 9 = 45°. 



where 



C"n..m, — 



2n + 1 (n — m)\ 



Air (n + m)\ 

The associated Legendre polynomials P™(cos#) have the following property: 

Vm ^ P™(±1) = 

^°(±i) = (±ir 



(4) 

(5a) 
(5b) 



Eq. [5] demonstrates that the field distribution along the z-axis, corresponding to cos# = 
±1 can be entirely described by taking the scalar spherical harmonics Y n $. Fig. [4ji shows the 
field reconstructed by imposing m = in Eq.Q with the same conditions of illumination 
as in Fig. [2|). Comparison of these two maps illustrates that the longitudinal in-phase 
resonances can be fully described when imposing m = and that this method can be 
considered as a relevant way to study the longitudinal modes of the structure. It must 
be stressed that these longitudinal contributions include the three directions of the field 
^n,o> ^n,o an d ~^ n ,o- The contribution of the electric field oriented along the chain axis 
is then given by the terms E^(r)Y n $ in Eq.(l). Fig. ij) shows the phase of the field 
reconstructed with the terms E^(r)Y n $. The in-phase oscillations of the longitudinal mode 
in the structure are clearly highlighted. Let us remark that the polarization moment inside 
the particles shows a clear dipolar-like behaviour. 

We now apply this formalism to study the phenomenon observed in Fig. [3j Fig. 
displays the longitudinal multipole contributions to the field at A = 595 nm and 6 = 45° 
obtained by reconstructing the field intensity map with (n,m) = (n, 0). It highlights that 
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FIG. 4: (a) Electric field intensity ||i?|| 2 /||i^ nc || 2 map in log scale of the longitudinal contributions 
reconstructed imposing (n, m) — (n, 0) in Eq. ([I]) and (b) phase of the field reconstructed with the 
terms E$(r)Y nfi in Eq. jl]) at A = 665 nm and 9 = 90°. 

the vanishing of the near field in cavity 1 associated with the high intensity enhancement 
in cavity 2, which are the striking features observed in Fig. [5]2, are essentially obtained. It 
shows that the intensity contrast observed is due to the excitation of longitudinal modes 
in the structure. More precisely, the phase of the longitudinal electric field displayed in 
Fig. [5^, obtained under the same conditions as Fig. shows that the mode L2 depicted 
in Table. |T] is predominantly excited. The oblique illumination has enabled the coupling 
of the incident light to an opposite-phase mode that is responsible for a vanishing of the 
field in cavity 1 while the field in cavity 2 is strongly enhanced. Let us emphasize that the 
electric field oscillations are principally dipolar in nature and that the longitudinal mode 
observed in Fig. [3] involves predominantly dipolar contributions. Unlike the case illustrated 
in Fig. |4]2, Fig. clearly shows a weak longitudinal induced polarization inside sphere 
1. The asymmetry of the electric field amplitudes inside the spheres in Fig. (5ja indicates 
not only that the L2 mode is excited, but that a combination of the three longitudinal 
modes (Table. [T]) is involved. This result shows that it is possible to nearly extinguish 
the longitudinal moment of a metallic nanoparticle by simply coupling this particle with a 
dimer of identical particles. It means that optical transparency [18J can be achieved with 
particles which do not support transverse modes, i.e. ellipsoidal or rectangular nanoparticles. 
Consequently, this study demonstrates that the dipolar moment of a plasmonic particle can 
be fully extinguished in ultra compact symmetric structures. 

This study demonstrates that antisymmetric modes can be excited in a linear trimer 
of identical particles by simply tuning the angle of incidence. Furthermore, we show that 
in a narrow frequency range, symmetric and antisymmetric mode-balancing extinguishes 
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FIG. 5: (a) Electric field intensity \\E\\ 2 /\\E inc \\ 2 map in log scale of the longitudinal contributions 
reconstructed imposing (n, m) = (n, 0) in Eq. ([!]) and (6) phase of the field reconstructed with the 

at A = 565 nm and 9 — 45°. 

the dipolar longitudinal moment of a metallic nanosphere coupled to a dimer composed of 
identical nanospheres. At this frequency, the total longitudinal moment of the trimer is 
weak and this phenomenon is not associated with an enhancement of the scattering cross- 
section of the trimer. This nearly zero longitudinal moment of sphere 1 leads to a near 
vanishing light intensity in cavity 1, and the high longitudinal moment of sphere 3 leads to 
strong enhancement of light intensity in cavity 2. Finally, we wish to emphasize that this 
phenomenon occurs in a simple fully sub- wavelength system by employing straighforward far- 
field control in the harmonic regime. Specifically, it enables the control of light localization 
between two identical nanosites separated by a distance of 55 nm, leading to a spatial 
resolution of A/ 10. 

The authors thank S. Bidault and R. C. McPhedran for stimulating discussions. 
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